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Abstract 

A preliminary study of time dependence of B^B^ oscillations using dilepton events is presented. The 
flavor of the B meson is determined by the charge sign of the lepton. To separate signal leptons from 
cascade and fake leptons we have used a method which combines several discriminating variables 
in a neural network. The time evolution of the oscillations is studied by reconstructing the time 
difference between the decays of the B mesons produced by the T(4S) decay. With an integrated 
luminosity of 7.7 fb~^ collected on resonance by BABMi at the PEP-II asymmetric B Factory, we 
measure the difference in mass of the neutral B eigenstates, Am^o, to be (0.507 it 0.015 it 0.022) x 
10^2 /is-^ 
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1 Introduction 



A precision measurement of the B^B^ oscillation frequency is of great importance since it is sensitive 
to the CKM matrix element \Vtd\ and, in combination with knowledge of the BgBg oscillation 
frequency, provides a stringent constraint on the Unitarity Triangle. 

The mass difference Am^o between the two mass eigenstates of the B^B^ system may be 
measured by comparing the rate for pairs of neutral B mesons to decay with the same b quark 
flavor with the rate to decay with the opposite flavor sign at the T(45') in the following time 
dependent asymmetry: 

N{B^B^){^t) - {N{B^B^){M) + N{B^W){M)) ^ 

N{BOB^){At) + {N{BOBO){At) + N{BOBO){At)) ^ bo ), {) 

where At is the difference between the two B meson decay times in the T{4:S) center of mass system. 
The simplest way to determine the b quark flavor of the decaying neutral B is to use leptons as 
tagging particles. By counting the number of "like" events + and "unlike" events 

a measurement of Am^o may be extracted through the asymmetry : 

, nA,n _ Njl^^n - iNil+,l+) + Nii-,i-)) 

The semileptonic (muon or electron) branching ratio of B mesons is about 20%. Therefore, 
the dilepton events useful for this analysis represent 4% of the T(4S) BB decays. In statistical 
terms, the dilepton tagging is more efficient than the semi-exclusive tagging performed at the 
ARGUS and the CLEO [Q] experiments. Moreover the new asymmetric B factories, like PEP- 
II, allow a time-dependent measurement, which is radically different from measurements of the 
time-integrated probability Xd performed at the previous e~^e~ colliders operating at the T(45), 
where Xd = x'^/ (2 • (1 + a^^)) and Xd = Am^o /T^o. Previous measurements of the time-dependence 
of B^B^ oscillations have been done by the LEP, SLD and CDF experiments Q . 

The present measurement is performed on events collected by the BABAR detector at the PEP-II 
asymmetric B Factory between January and June 2000. The corresponding integrated luminosity 
is 7.7 fb"^ taken on the T{AS) resonance and 1.2 fb^^ taken 40MeV below resonance. The BABAR 
detector and its performance are described elsewhere [^. The event selection and particle identifi- 
cation criteria are described in section The selection of signal events and a study of the fraction 



of events with the wrong flavor tagging (mistag) are detailed in Section 2A. The method to deter 
mine the time-separation of the two B semileptonic decays is explained in Section ^. Section 4.2 
shows the details of the fit on data and the result of the Am^o measurement. A list of cross-checks 
of the result is in Section |5|, while the evaluation of systematic uncertainties is reported in Section 

H 



2 Selection of dilepton events 

In this study of the oscillation frequency Am^o , the flavor of the B meson at decay is determined 
by the sign of leptons produced in semileptonic B decays. To reduce the mistag rate, an attempt 
is made to suppress cascade leptons (produced in 6 — > c ^ £ transitions). 
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2.1 Lepton identification 

Electron and muon candidates are required to pass the very tight selection criteria fully described in 
Q. Electrons are selected by specific requirements on the ratio of the energy deposited in the Elec- 
tromagnetic Calorimeter (EMC) and the momentum measured in the Drift Chamber (DCH), on 
the lateral shape of the energy deposition in the calorimeter, and on the specific ionization density 
measured in the DCH. Muons are identified by the use of the energy released in the calorimeter, 
as well as the strip multiplicity, track continuity and penetration depth in the Instrumented Flux 
Return. The performance of the very tight selection criteria are estimated on data control samples, 
as a function of the particle momentum as well as the polar and azimuthal angles. The electron 
and muon selection efficiencies are about 92% and 75%, respectively, with pion misidentification 
probabilities around 0.3% and 3%, respectively. Lepton candidates consistent with the kaon hy- 
pothesis as measured in the Detector of Internally Refiected Cherenkov light (DIRC) are rejected. 
More than 60% of the kaon contamination in the muon sample is rejected with negligible effect on 
lepton identification efficiency. 

2.2 Background rejection 

Non BB events are suppressed by requiring the Fox- Wolfram ratio of second to zeroth order mo- 
ments to be less than 0.4. 

The residual contamination from radiative Bhabha and two-photon events is reduced by requir- 
ing the event squared invariant mass to be greater than 20 (GeV/c^)^, the event aplanarity to be 
greater than 0.01, and the number of charged tracks to be greater than 4. 

Electrons from gamma conversions are identified (see @) and rejected with a negligible loss of 
efficiency for signal events. Leptons from J/il^ decays are identified by pairing them with the other 
oppositely-charged candidates of the same lepton species, selected with looser criteria. We reject 
the whole event if any combination has an invariant mass within 40 MeV/ of the J/i/^ mass. 

2.3 Track quality requirements 

We finally apply selection criteria on the quality of the tracks, in order to improve the Az recon- 
struction, where Az is the difference between the decay points of the two B mesons along the beam 
direction. Any lepton candidate must have a distance of closest approach to the nominal beam 
position in the transverse plane, do, less than 1 cm, and a distance of closest approach along the 
beam direction, \zq\, less than 6cm, at least 20 hits in the DCH, at least 4 z-coordinate hits in 
the Silicon Vertex Tracker, a momentum range in the center of mass system between 700MeV/c 
and 2.5GeV/c, a momentum range in the laboratory system between 500MeV/c and 5GeV/c, and 
a polar angle in the range between 0.5 and 2.6 radians. We also require the total error on Az, 
computed on an event-by-event basis to be less than 175 /im. When estimating the event-by-event 
error, it should be noticed that, due to non-zero flight length of the B mesons in the transverse 
plane, the two leptons do not actually originate from the same point in that plane. The total error 
is therefore the quadratic sum of the tracking error and of this additional uncertainty. As reported 
in Section ^, the non-negligible effect of the track quality requirements on signal efficiency yields 
only a small degradation of the resulting statistical uncertainty in Am . 
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2.4 Selection of the direct dileptons 

The discrimination between direct and cascade leptons is based on a neural network which combines 
five discriminating variables, all calculated in the T(4S') center of mass system: 

• the momenta of the two leptons with highest momenta, and 

• the total visible energy, Etot, and the missing momentum, Pmiss, of the event; 

• the opening angle between the leptons, 012- 

2000 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 
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12 3 

Figure 1: Distributions of the discriminating variables (a) p*, (b) P2, (c) Etot, (d) Pmiss, and (e) 
612, for data (points) and Monte Carlo (histograms). The contributions from direct-direct pairs, 
direct-cascade pairs, and pairs with one or more fake leptons, are shown for the Monte Carlo 
simulation. 

The distributions of these variables are shown in Figure |l], for data and Monte Carlo simulation. 
The first two variables, pi and P2, are very powerful in discriminating between direct and cascade 
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leptons. The last variable, 9 12, efficiently removes direct-cascade lepton pairs coming from the same 
B and further rejects gamma conversions. Some additional discriminating power is also provided 
by the other two variables. The chosen neural network architecture (5:5:2) is composed of 3 layers, 
with 2 outputs in the last layer (one for each lepton). The network is trained with 40,000 dileptons 
from generic 5° and B^, and the outputs are chosen to be 1 and for direct and cascade leptons 
respectively. Figure |2| shows good agreement between data and simulation for the neural network 




Figure 2: Neural network outputs distributions for the (a) highest and (b) second- highest momen- 
tum leptons, for data and Monte Carlo simulation. The various Monte Carlo contributions are 
shown separately. 

outputs of the two leptons. We require both outputs to be greater than 0.8. 

The combined effect of the above cuts gives, from simulated events, signal purity and efficiency 
of 78% and 9%, respectively. The remaining background consists of 12% direct-cascade events (8% 
with the wrong tagging), 5% BB events with one or more fake leptons, 2% BB events with one 
or more non-prompt leptons, a negligible contribution from cascade-cascade events, and 3% from 
continuum events. The latter was determined in data by rescaling the number of off-resonance 
events that pass the selection with the ratio of on- and off-resonance luminosities. The total 
number of selected on-resonance events is 36631 (10742 electron pairs, 7836 muon pairs, and 18053 
electron- muon pairs). 

3 Determination of At 

A determination of the z coordinate of the B decay vertex using only the lepton track can be 
obtained, to first approximation, by taking the z of the point of closest approach between the track 
and the beam spot in the transverse plane. This estimator is a fairly good way to determine the 
z position of the B^ decays vertices since the selected direct leptons have rather high momenta. 
However, it is possible to use the two lepton tracks and a beam spot constraint in a simple 
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vertex fit to obtain a better estimate of the primary vertex of the event in the transverse plane, and 
to compute the points of closest approach of the two tracks to this new point. The corresponding z 
coordinates represent a better approximation of the z coordinate of the B decay vertices, and the 
corresponding Az resolution function has much reduced tails. For this reason we adopt the latter 
method to compute the Az. 

Further studies show that a requirement on the total error to be less than 175 /xm reduces 
the tails of the Az resolution function by a factor four and reduces the signal efficiency by 30%. 
However, due to the improved resolution, the total statistical uncertainty on Am^o is degraded only 
by 3% despite the loss of efficiency. A two-Gaussian fit to the resulting Az resolution function from 
simulated dilepton events gives = 87 iiva. and = 195 /xm for the narrow and wide Gaussian, 
respectively, and 76% of the events in the narrow Gaussian. 

The time difference between the two B decay times is defined as At = lS.z/{< /Sj > c), with 
< /37 >= 0.554. This approximation neglects the B meson motion in the T(45) rest frame. In 
this inclusive approach it is not possible to determine the exact boost. Therefore, the effect of this 
shift was studied with Monte Carlo by comparing the fitted value of Am^o with the true At and 
with Az/{< (3'j > c). This study shows that the effect is negligible compared to the current level 
of accuracy of this analysis. 

4 Fitting procedure 

4.1 Time dependence of the fraction of mistagged events 

Even after a cut on the neural net output, a non- negligible fraction of events are mistagged (i.e. 
a true B^B^ pair is tagged as a B^B^ or B^B^ pair and vice versa for B^B^ or B^B^ events). 
The fraction of mistagged events is directly determined in the fit. However we have to take into 
account that the time dependence of cascade leptons from the same i?, or from the other are 
different. In the case of a cascade from a same B (Fig. ^(a)), we observe a peak at low Az, due 
the fiight length of the charm hadron, which is fitted by an exponential decay. For the cascade 
leptons from the other (circles in Fig. |3|(b)), the fraction of mistagged events as a function of 
the true Az between the two leptons shows a linear dependence which comes from the fact that the 
Az measured between the two leptons contains the additional flight length of the charm hadron. 
The same distribution, determined using the z distance of the true B vertices, is flat to flrst order 
(squares in Figure]^. Actually, the linear dependence on Az can be explained by considering the 
time distribution of the cascade lepton. Assuming that the fiight length of the charmed hadron is 
small compared to that of the B, the time distribution of the cascade lepton from the other B can 
be approximated by: 

^e-t/(rs+<r.» ^ ^(1 + ±^ . t)e-*/-fl, 

where tb and < Tc > are the B meson lifetime and an average lifetime of the D mesons, respectively. 
In the final fit procedure, the linear dependence is taken into account by a free parameter; the time 
shape of the cascade lepton from a same B is determined from the Monte Carlo simulation. 

4.2 Measurement of Athbo 

The value and statistical error for Am^o are extracted with a minimization fit to the dilepton 
asymmetry (see Eq. |2|). The fit function, Afit{At), takes into account the various time distributions 
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Figure 3: (a) Distribution of Az between direct and cascade leptons that come from the same B. 
(b) Fraction of mistagged events for cascade leptons coming from the other i3, as a function of the 
true Az between the two leptons (circles) and the true Az between the B mesons (squares). 



of the dilepton signal (/""'™*^(At), /™*^(At)), the cascade lepton and the non-BB backgrounds 
(/o^L(At), /f,L(At)): 



, .... if' - f') ® fresojAt) + (/OL(At) - &,m) 
'"'^ ^ U^' + P') ® fresoiAt) + (/«L(A^) + fofneri^t)) ' 

where stands for the convolution product with the resolution function freso{At) (see Section ^), 
and the /-functions are expressed in terms of the various signal and background contributions as 

f^{At) = f~{At)-ii-{i + fM + &a,i^t)-fc-m 

+n^{At).{l^, + aAt) 
f^{At) = /'"-(At).(l-(l + /e)%) + /— (At).(7?o + aAt), 



The signal contributions for unmixed and mixed events are given respectively by 

f~^{M) = — i— -[rV^"l^*l(l + cos(AmBoAt)) + 2i?-r+-e-^''l^*l], 
2(1 + R) 

^""^^^^ = ^^^^[rVr"l^*l(l-cos(AmBoAt))], 

where R is proportional to the ratio {h\f-\ )/(6o/oo) (&+ and 6o are respectively the semileptonic 

branching ratio of charged and neutral B, and /-| //oo is the production ratio of charged and 

neutral B pairs at the T[AS)). 
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The time distribution of direct-cascade events where both leptons originate from the same B 
is represented by fmfstaqi^^) = [< T'^ > e~"^'"''^l^*l], with < c/V^ >= 60 ;um as determined from 
simulated events (Fig. pi). The difference in the fraction of direct-cascade events between the 
cascade lepton from the same B and from the other B (due, for instance, to the cut on the angle 
between the 2 leptons) is estimated by the parameter /c, determined to be 0.6 from simulated 
events (Fig. ^). 

The time dependence observed for the mistag fraction of direct-cascade events where the cascade 
lepton comes from the other B, as discussed in Section 41, is parametrized by a constant term, rjQ, 
and a slope, a. The same functional dependences as for signal events, fu'^"^^^(^/\f'j and /'"*^'(At), 
are used. 

The time distributions of the non-Si? background, /^^.^.(At) and /^^g^(At), and their absolute 
normalizations, are obtained from off-resonance data. 

Four parameters, Am^o, r/o, R, and a, are fitted directly to the observed asymmetry. The 
lifetimes of the charged and neutral B, and F", are fixed to their world average values |^. 

The off-resonance data is used to measure the fraction of the non-BB background to be (0.7 ± 
0.1)% for same-sign dileptons and (2.2 it 0.3)% for opposite-sign dileptons, respectively. A fit to 
the time distribution of these events yields an effective lifetime equal to 130 /im and 135 ^um for 
same-sign and opposite-sign dileptons, respectively. 

The fit to the measured asymmetry Afui^t) shown in Fig. ^ and obtained with an integrated 
luminosity of 7.73/6"^, yields the following values: Am^o = (0.507 ± 0.015) x lO^^Tis"^, % = 
0.109 ± 0.004, R = 1.34 ± 0.11 and a = (-1.7 ± 3.3) x 10-^ with a of 20.8 for 21 degrees of 
freedom. 




Figure 4: Distribution of the measured asymmetry 74of,s(|At|) between unlike-sign events 
and like-sign events + for (a) the inclusive dilepton sample and (b) the dilepton 

sample enriched with soft pions, which is discussed in Section The curve represents the result of 
the fit. 
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5 Cross-checks and stability of the Am^o measurement 



5.1 Enrichment of the neutral B with a soft pion 

In the inclusive approach proposed in this analysis, the final dilepton sample contains both charged 
and neutral mesons B in almost equal proportions. Therefore, the observed oscillation amplitude 
is reduced by the presence of the non-oscillating charged B. To enrich the B^ fraction, the direct 
lepton can be correlated with the soft pion produced by a D*~^ decay. Charged B mesons can only 
produce a direct lepton and a charged D* through the D** decay or through the non-resonant 4 
body decay B~ 0*^71^ i~V. The branching fractions of these modes are not perfectly measured, 
but they should represent roughly 10-20% of the semileptonic decays. 

The identification of an event with a soft pion is based on a method proposed by the CLEO 
Collaboration only tracks with momentum less than 190MeV/c in the center-of-mass system 
are considered. The direction of motion of the D* is very close to that of the soft pion (the and 
the soft pion are produced almost at rest in the D* system) and the energy E"^* of the D* in the 
T{4:S) system is approximated by using the energy of the soft pion E* in the T(4S') system and 
the energy of the soft pion E!f^* in the D* system: E*j~), ~ {E*/E^*) ■ Md*. With the four- vector 
of the lepton and the D*, one can compute the missing mass squared of the neutrino. In the 
analysis, events are kept if |M^| < 1.0( GeV/c^)^. 

The fit of this sub-sample gives Am^o = (0.518 ± 0.017) x 10^^ hs^^, in good agreement with 
the value obtained with the dilepton sample. Even though the fraction of events with the additional 
soft pion represents only 16.5% of the total dilepton sample, the statistical errors are comparable. 
While, for the moment, this preliminary method constitutes an excellent cross-check, it may later 
become an alternative approach in its own right. 

5.2 Stability studies 

We have investigated the stability of the fit results against various changes in selection criteria. The 
fit was performed in several ranges of azimuthal angles, as well as for a range of values for the cut 
on the neural network outputs, (0.6 0.9), on the total error of the Az (150 /im ^ 300 /im) and 
for a range in Az. Subsamples composed of only fifi, ee and were also considered. In all cases, 
variations in Am^o were found to be small or consistent with the nominal value within statistical 
errors. 

6 Systematic uncertainties 

In this analysis the fraction of mistagged events % is directly extracted from the fit of the asymmetry 
but a time dependence of this component, as well as the fraction of misidentified leptons, may induce 
a bias in the Am^o determination. These effects are corrected by using the time distribution of the 
mistagged events determined from the Monte Carlo, and by fitting a slope to the mistag fraction 
time dependence. The systematic error is determined by assuming that the Monte Carlo corrections 
are known at the 30% level. 

A conservative estimate of the uncertainty due to fake leptons is taken to be the difference 
between the results of the fit to Monte Carlo with perfect and simulated particle identification. 

Another important source of systematic errors comes from the determination of the resolution 
function, which is taken from simulated events. To estimate the uncertainty involved in this proce- 
dure, we have compared the Az resolution between data and Monte Carlo using J/'i/' events, where 
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the leptons are known to come from the same vertex. From this comparison, shown in Figure ^, we 
estimate an uncertainty on the width of the narrow and wide Gaussians of the resolution function 
of 5% and 10%, respectively. 




Az (i^m) Az (i^m) 

Figure 5: Az distribution for lepton pairs in the J/iIj mass window in (a) data and (b) Monte 
Carlo events. The distributions are fitted with a sum of two Gaussians. The resolutions for the 
narrow and wide Gaussians are 101 /im and 205 ^m in data, 102 /im and 184 //m in Monte Carlo, 
respectively. 

The effect of a charge asymmetry in the identification of the lepton (e+ ^ e~) or a mistag asym- 
metry 1]'^ 7^ rj^ on the Am^o measurement is negligible since the effects cancel in the asymmetry. 
However, the mistag probability r] may be different for the charged and neutral B. The impact of 
such an effect on the Am^o measurement is negligible because the bias is fully absorbed by the 
parameter i?, which implies that the fitted value of this ratio need not necessarily be unity. 

The list of systematic effects is summarized in Table ^. The sum of the different contributions 
gives a total systematic uncertainty of 0.022 x 10^^ fis^^. 

7 Conclusions 

We present a preliminary study of the B^B^ oscillation frequency with an inclusive sample of 
dilepton events corresponding to a total luminosity of 7.73 fb~^ collected by the BABAR. experiment. 
We obtain Am^o = (0.507 ± 0.015 ± 0.022) x 10^^ h s~^. The accuracy is already comparable with 
the current world average. 



16 



Table 1: Summaxy of the contributions to the systematic uncertainty in Am^o. 





(T(Amgo) 


or 

Source or systematic uncertainty 


(-1 rii2 t 

(1\J lis ) 


TV T T ") T ") T 1 -1 

JN on- is is background 


0.005 


Mis-identihcation 


0.011 


Time-dependence or the cascade events 


0.009 


Correction or the boost approximation 


0.001 


y-motion of the beam spot 20 /xm) 


n nni 

U.UUi 


Az resolution function 


0.009 


Tails of the Az resolution function 


0.004 


Time-dependence of the resolution function 


0.006 


Sensitivity to r+ (PDG 98 ±la) 


0.007 


Sensitivity to (PDG 98 ±la) 


0.007 


Total 


0.022 
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